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The Keplerian motion of accretion disks in Active Galactic Nuclei 
(AGN) is usually believed to be generated by a heavy central mass. We 
investigate accreting disk systems with polytropic gas in Keplerian rotation 
and obtain a phenomenological formula that relates the Keplerian angular 
frequency to the ratio of disk and central masses. Central mass approaches 
the Keplerian value, if the inner boundary of a disk is close to the minimal 
stable orbit of a black hole. These results are applied to NGC 4258, the 
unique AGN with a finely measured Keplerian rotation curve of the central 
disk, with the conclusion that its rotation curve is, in fact, determined by 
the central black hole. The mass of the accretion disk exceeds IOOMq. 

PACS numbers: 95.30.Lz, 95.30.Sf, 98.62.Js, 98.62.Mw 

In the astro-physical literature, the velocity of rotation of disks in Kep- 
lerian motion is typically associated with a heavy central mass M c and a 
light disk mass The discovery of the massive black hole in the masing 
nucleus of NGC 4258 exemplifies this interpretation [1] . On the other hand, 
it is known that even the selfgravity of heavy disks can be compatible with 
the Keplerian rotation [2]. Hure et al. [3] investigated thin dust disks and 
pointed out that the rotational Keplerian velocity gives only information on 
the enclosed mass of the system, and it is not possible to separate masses 
of the two components — of a disk and a central mass. 

We show in this paper that the Keplerian rotation curve, with the fre- 
quency id, does depend on the two masses but they do not contribute addi- 
tively. A phenomenological formula expresses M c G/(r 5 uj 2 ) (r and G are the 
cylindrical radius and the gravitational constant, respectively) as a function 
of x = Md/M c and of the ratio y = r- m /r out , where r- m and r out are the 
innermost and outermost disk radii, respectively. If y is small enough then 
the central mass is reasonably well approximated by the Keplerian value 
r 3 oj 2 /G. 
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Consider a stationary, axially symmetric, selfgravitating disk of gas, ro- 
tating around a central point mass M c . The gravitational potential of the 
system can be written as the superposition $ = —GM c /\x\ + $ g . The po- 
tential $ g is due to the gravity of the disk and satisfies the Poisson equation 

A$ g = AirGp, (1) 

where p denotes the mass density of the gas. In the following, we use the 
cylindrical coordinates (r, (ft, z). The velocity of the gas reads U = oo{r, z)d ( j ) . 
For a barotropic equation of state the frequency oj depends only on the 
cylindrical radius r [I]. We assume polytropic equations of state p = Kp r , 
where p is the gas pressure, and K and V are constant. The Euler equations 

Vp + p(U • V)U + /?V$ = 

can be integrated, yielding 

h + $ c + $ = C (2) 

in the region where p ^ 0. Here, h denotes the specific enthalpy of the fluid 
dh = dp/ p, and 



/r 
dr'r'oj 2 (r') 



is the centrifugal potential. The structure of the disk can be obtained from 
Eqs. ([I]) and Q provided that the equation of state and the rotation law 
oj = oj(r) are known. 

Let the Keplerian rotation law be w = oj^/r^/ 2 . It follows from Eq. Q 
that 

= GM C + nnTWt ~ * g(nn)) , (3) 

?"out '"in 

where $ g (r in ) = $ g (r = r in , z = 0) and ^ g (r out ) = $ g (r = r out , 2 = 0). We 
fix r out and the maximal density /? max - Then, we find that the limit of the 
right-hand side of Eq. ([3]) for r- m — > is the Keplerian angular frequency 
Wq = GM C . To show that, note that <3? 9 is a bounded function on M 3 . It is 
also possible to prove that oj 2 — > GM C as r; n — > r ou t- 

In the numerical analysis we apply the Self-Consistent Field (SCF) 
scheme [5J El [7] . Numerical solutions converge quite well for Keplerian ro- 
tation, although the drawback of the method is the need to employ a large 
number / of Legendre polynomials for thin disks. We find accurate solu- 
tions with / = 400 and large grids (up to 5000 x 5000). Such resolutions 
are permitted because of the simplicity of the SCF method. Our numerical 
solutions have been tested by the new version of the virial theorem that in- 
cludes a point mass [8]. It should be pointed out, however, that even I = 400 
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Fig. 1. The dependence of u>q/(GM c ) on the mass ratio x. The plot is obtained for 
the adiabatic index T = 5/3 and Hn/^out = 10 -4 . The dashed line depicts analytic 
fit Q. Crosses correspond to numerical solutions obtained for different values of 
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Fig. 2. Mass density profile within a disk model obtained for T = 5/3. The units 
are g cm~ 3 . 
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does not suffice to resolve structures in the gravitational potential that are 
thinner than thinnest disks considered in this paper. We performed tests on 
thicker disks, demonstrating that important quantities are insensitive to the 
addition of a higher number of expansion polynomials, which could model 
fine structures of the gravitational potential within the disk. This conclu- 
sion is also supported by the heuristic argument provided by the analysis 
of infinitely thin disks — they give a regular gravitational potential, i.e., 
without any structure comparable with the disk thickness. 

We found a few thousands of numerical solutions for polytropic disks, 
mainly with two polytropic indices T = 5/3 (which can be thought to rep- 
resent monoatomic gas) and T = 4/3 (which might well be associated with 
vapor water maser disks) and for a few dozens of values of the interior ra- 
dius corresponding to y G (l0 _4 ,0.5). Calculations have been done also for 
a range of polytropic indices in between 4/3 and 5/3, and the internal radius 
corresponding to y = 10 -4 . 

Below we shall report results concerning the case T = 5/3, but the 
essential features of solutions are similar in other cases. The obtained results 
can be quite well approximated in the quadrant < x < 1, 10~ 5 < y < 0.55 
by the following phenomenological formula 

<"° -! + „(,) 2M±£, (4) 



GM C JKyj h(y) + x' 

where x = M^/M c and /, g and h are defined as follows 

f(y) = 1.322y - 3.719y 2 + 13.47y 3 - 44.48y 4 

+116.9y 5 - 206.2y 6 + 208.9y 7 - 90.60y 8 , 

g(y) = 0.4778 + 0.9495y - 31.67y 2 + 258.4y 3 

-1431y 4 + 5169y 5 - lllOOy 6 + 12760y 7 - 6034y 8 , 

h(y) = 0.7423 + 0.6313y - 43.57y 2 + 375y 3 

-2059?/ 4 + 7281y 5 - 15390y 6 + 17520y 7 - 8242y 8 . 

Functions /, g and h are nonnegative and smaller than 1 for arguments 
< y < 1. 

It is clear from these formulae that the central mass is always smaller 
than Lu 2 r 5 /G. At the same time, it is easy to see from Q that if y is 
sufficiently small and Md < M c , then the rotation curve of the gaseous disk 
is almost exclusively influenced by the central mass. 

The number of known AGN's with masers is increasing in recent years 

[aEgimEananiiEgnB] ( see also for a review us )> but there there is 

only one case of well resolved Keplerian rotation curve — in the AGN of 
NGC 4258. Observations of supermasers near the center of the galaxy NGC 
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4258 [TJ [18] provide evidence for a rotating disk of gas surrounding a black 
hole with a mass M c = (3.9 ± 0.3) x 10 7 Mq. The motion is Keplerian to 
high precision, of about a part in a hundred [U \19\ \20\ I2T] , and this mass 
estimate assumes that velocities are due solely to the gravitational force of 
the central body. A similar result has been obtained by Siopis et al. |22j . 
who determined the mass of the black hole by constructing axisymmetric 
dynamical models of the galaxy NGC 4258. Their best mass estimate yields 
M c = (3.3 ±0.2) x 1O 7 M . 

The masers are observed in the annular region that extends from 0.13 pc 
to 0.26 pc from the center, and the region with detected masers is thin - 
the maximal relative height is circa 1 /800 (these distances can be expressed 
in terms of the Schwarzschild radius: r g = GM c /c 2 = 5.8 x 10 12 cm and 
lpc = 5.4 x 10 5 r g ). There exists evidence that the disk fills not only the 
annular region in which masers are detected, but that it extends outwards 
from a vicinity of the black hole horizon. Reynolds et al. [23] investigated 
the weakly radiating AGN in NGC 4258 using data from Suzaku, XMM- 
Newton, and the Swift/Burst Alert Telescope survey. They constrained its 
luminosity region to between 10 and 4 x 10 4 Schwarzschild radii of the black 
hole. In the numerical calculations we assume that the disk extends from 
r; n = 2.6 x 10~ 5 pc to r ou t = 0.26 pc, and its relative height is equal to 
1/800. 

It is clear from (jij) — see also Fig. 1 — that if M d < M c then u%/(GM c ) « 
1 with the deviation from unity much smaller than 0.001. The disk's mass, 
obtained numerically for the polytropic index T = 5/3 and assumed geo- 
metrical parameters of the disk, reads M^ = 2.62 x 10 _6 M C ~ 102Mq. 

The results for other polytropes, 5/3 > T > 4/3, suggest that disk's 
mass depends rather weakly on the equation of state and increases (up to 
Md = 190Mq) with the decrease of the polytropic index. The disk's size 
can be larger than that detected in observations, which, in turn, would yield 
a larger mass. Therefore M<j = IOOMq can be regarded as a conservative 
lower bound on the mass of the gaseous disk in the AGN of NGC 4258. 
A sample distribution of the mass density within the disk, obtained for 
r = 5/3, is shown on Fig. 2. 

The above results are consistent with those rough mass estimates that 
base on such values of the baryonic number density p n within the disk that 
favor the existence of water masers. If one assumes that p n is around 10 10 
per cm , then the mass is around 1O 4 M [HQS]. Herrnstein et al. [21] note 
that the allowed density is 10 8 /cm 3 < p n < 10 10 /cm 3 ; the lower bound gives 
the disk mass ~ IOOMq. 

This analysis is entirely Newtonian, and since it applies to systems with 
black holes, it should be repeated in the general-relativistic context. We do 
not expect essential changes concerning those systems, where the inner disk 
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boundary is well separated from a black hole horizon, r- in 3> 6M C . 
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